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Abstract 

Experimental projects using spherical antennas to detect gravitational waves are nowdays a 
concrete reality. The main purpose of this paper is to give a possible way of interpreting output 
data from such a system. Responses of the five fundamental quadrupole modes and of the six 
resonators in TIGA collocations are shown as a function of the incoming direction of the incident 
wave. Then, for a source lying in the galactic plane, sidereal time and galactic longitude dependence 
is given. Thus, once a candidate source of gravitational waves is considered, we can exactly predict 
the resonators' response as a function of time. 



* Electronic address: |alice.gasparini@physics.unige.ch| 



1 



I. INTRODUCTION 



Theoretical interest in spherical gravitational wave (GW) antennas dates back to the 
1970s. Spherical antennas have a greater cross-section than a bar detectors of similar dimen- 
sions. More importantly, they have both omnidirectional and omnipolarization sensitivity, 
and also the potential to detect the direction of wave provenance. 

Interest in experimental research into resonant spheres has increased over the past 15 
years, and today spherical antennas are recognized to be the new generation of gravitational 
resonant detectors, to complement existing cylindrical antennas. Two experiments are under 
way: MiniGRAIL in Leiden (Holland) |l| and The Graviton Project in Sau Paulo (Brasil) 

One of the main problems today in detecting gravitational waves with resonant spheres 
is to find the "best" location for amplifiers on the detector surface in order to optimize 
sensitivity. To date, no theoretical basis has been developed to provide a definitive solution 
to this problem. Nevertheless many suggestions have recently been proposed, including 
the uncoupled transducer configuration by Zhou and Michelson in 1994 [3], and the PHC 
configuration by Lobo and Serrano in 1996 B; both are five- resonator configurations, the 
first including radial and tangential motion of the resonators, the second only radial motion. 
In addition, the TIGA configuration, proposed by Merkowitz and Johnson 5] (1995), consists 
of six transducers moving radially, and it enjoys a peculiar symmetry which greatly simplifies 
equations. 

The general theoretical aspects of an ideal simple sphere interacting only with gravita- 
tional radiation have already been treated in the literature (see for example |7] or Ig]). In 
the first part of this paper we will work out briefly some of them and obtain the energies 
stored in the five fundamental sphere quadrupole modes, as well as their dependence on the 
incident wave direction and polarization. 

Taking into account of earth's rotation and the detector's location, we will translate the 
direction dependence into galactic longitude and sidereal time dependence for a source lying 
in the galactic plane. 

We will review the general problem of the spherical detector with radial moving resonators 
on its surface, again in the case of high signal-to- noise ratio (SNR). We then focus our at- 
tention on the TIGA configuration, finding the oscillation amplitudes for the six transducers 
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function of the incident GW direction. 
With the same method used in the case of the simple sphere without amphfiers, we will 
find the resonators' amplitude dependence on sidereal time and on galactic longitude for a 
source lying in the galactic plane which emits randomly polarized GW radiation. 

The paper is organized as follows. In section 2 we briefly sketch some general aspects of 
the physics of GW resonant detectors without any assumption on the shape of the antenna. 
In section 3 we focus our attention on the spherical detector, and we present the solutions to 
the equations of motion followed by the five modes angular and sidereal time sensitivity. In 
section 4 results of the previous sections are generalized to the case of the spherical detector 
coupled to a set of radial transducers. Section 5 contains the summary and conclusions. 



II. GENERIC RESONANT DETECTORS 



A. The Equations of Motion 

We will work with a solid whose density is p(r), where r = (ri, r2, r^) = {x, y, z) is the 
position of an infinitesimal mass element relative to the center of the solid and we call the 
Lame coefficients of the material A and /i. 

We suppose that a gravitational wave hits our detector at the time t = 0. If we call 
^(r, t) the small displacement from the equilibrium position at the point r and time t, the 
equation of motion turns out to be |9| 

P^^^^ - PVMr, t)-{\ + p) V(V ■ u{r, t)) = Fir, t), (1) 
where F^(f,t) = pRlf,QX^{r, t) is the j^^ component of the tidal force density attributable to 



the GW 



id- 



As is well known for a resonant detector the spatial dependence of the Riemann tensor 
can be neglected. Moreover, as all our summation indices are spatial and we reasonably 
assume the background metric to be flat on the earth's surface ( g^u = rj^y + /i^j, ), we can 
write them all as lowered indices. Thus, in the first order of h, the Riemann tensor can be 
written -Rjojo = \hij{t), which is a rank 3 traceless symmetric tensor. It represents then a 
particle of spin 2, the graviton. 

Such a tensor can be decomposed into a base of 5 matrices. In order to do that, we follow 
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the procedure used by Lobo working with a basis of real matrices instead of imaginary 
matrices, as is usually done: 
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Noting / the identity matrix, the above matrices have the following properties: 
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in 



5^p, /,,M- = 0, 



with a, P = 0, Ic, Is, 2c, 2s. 

Then we can rewrite the j component of the force density: 
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Repeated indices are summed. Eq.(jll) introduces the functions 

/;(f) = pMf^Tk. 

Direct computation gives for the /°: 
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We now want to write down the explicit expressions for the functions h"{t) given a GW 
propagating along a generic direction with respect to the detector coordinates system. 
In the frame where the GW propagates along the z axis we can pass into the TT gauge in 
order to rewrite the Riemann tensor as 
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We call 9 the angle from the z axis and the angle between x and the projection into the 
X, y plane, {x, y, z) is the detector frame basis. For a generic arrival direction given by the 
angle {6, 0), we have to perform a rotation in order to obtain the Riemann in the antenna 
frame: 
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We now are able to compute the components h°'{t) which are 



(7) 



(8) 



sm\9)h+{t) 



[Att / 

W — (^- cos(6') sin(6') cos(0)/i+(t) + sin(6') sin(0)/ix {t 
(cos{9) sm{9) sm{(p)h+{t) + sin(6') cos(0)/ix (t 



^ (^(l + cos'^{9)) (cos\(f)) -l]h+{t)-2 cos(^) sin(0) cos(0)/ix (t) 
15 




((l + cos2(^)) cos(0)sin(0)/i+(t) + cos(^) (cos^{(p) - sin^icj))) /ix(t)) 
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Note that we have defined the + and x polarization with respect to the detector frame. 
In the case where the source has an intrinsic polarization e+,ex, we have to make the 
substitution 

h+ = cos(2'?/;)e+ + sin(2?/')ex 

hx = — sin(2?/')e+ + cos(2-?/')ex (10) 

where i/j describes the polarization angle. This becomes useful if we want to take the average 
of the polarization. 

B. The general solution 



The solution of Eq. (Q) can be expressed formally by means of the Green function integral 



getting: 



^ir,t) = T.T.^NYNQ%mN{^ (11) 

a N 

where 



= / $^(f)-r(r)rfV, 

J Solid 



Q^{t) = f h"{t') sin uJN{t-t')dt' (12) 
Jo 

and $Ar(r) are the normalized (i.e. they verifies /soUd (^) ' ^N'i^pi^d^r = M5]^^m') 
eigenfunctions solutions of the corresponding homogeneous equation 

pu;^$7v + /iV^$jv + (A + ;u)V(V-$7v) =0 (13) 

with the appropriate boundary conditions. indicates a collective quantum number denot- 
ing the quantum state. 

Putting the above result in a different form, i.e. writing the solution as 

M(f,t)=^i?^(t)$^(f), (14) 

N 

it is obvious, as expected from [l^ that the quantity = J2a ^N^fN^Ni^) satisfies equation 



Q JSolid 

= / M-^p{r)R,oko{t)rk^NAr). (15) 

J Solid 

Each mode is then formally equivalent to a one-dimensional harmonic oscillator with 
frequency u^, driven by a force per unit mass = J2afN^°'i'^)- -^^^ such a system,^ the 
expression for the energy per mass unit adsorbed from the driving force is given by [12]: 

1 f+OO 

E, = -\ FM{t)e-''^^'dt\\ (16) 

Z J — oo 

Before moving on to the case of the spherical detector, it would be intersting to use this 
general method on a cylindrical detector. While we omit this exercise here, it can easily 
be verified that this general method reproduces the results for the eigenfunctions and the 
energy stored in oscillations in the well-known case of the one- dimensional homogeneous bar 

Q. 

III. THE SPHERICAL DETECTOR 
A. The normal modes 

The normal eigenfunctions of Eq. pH|l for the sphere can be found in the literature (see, 
for example Q). Note that they are usually given in terms of imaginary spherical harmonics. 
Here we will put all results in terms of real spherical harmonics, but the shape of the solutions 
is the same. 

Let us briefly recall the main points. There are two kinds of solutions with boundary 
conditions crijUj = at r = R, with aij the stress tensor, conventionally defined as aij = 
^kkSij + ^fxuij, with Uij = ^{uij + Uj^i). These solutions are called the toroidal ones and the 
spheroidal ones: 

$f,„(r) = Ani{r)Yi^{e, <P)h - Bni{r)in x LYi^{e, </>) (17) 

where each mode n, /, a corresponds to the generic in Eq. ()13|) . n is a positive integer which 
represents the energy level for a fixed angular momentum I and Yia are the real spherical 
harmonics with kinetic momentum / obtained by the imaginary ones Yi^±rn in the following 
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way: 
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with m = 0, So defined, all the spherical harmonics ()18j) are real quantities. For / = 2 
this gives: 
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Also, given a generic unity vector in spherical coordinates, n 
we can see that 

Y2,^i6,<p) = Mt,n,nj 

where are the matrices defined in Eq.Q. 

Ani(r), Bniir) and Tniir) are scalar functions of r: 
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where the ji are spherical Bessel functions and 

ax X ax- 



Psix) = ^P2ix) + 



Poix). (22) 



The k and the q are related to the quantized u of Eq. (fT!?|) by 

kli = ^ qli = ^^. (23) 

B. Interaction with a gravitational wave 

We start by computing f^^^, of Eq.lfT^ for our spherical detector for both toroidal and 
spheroidal eigenfunctions: 



J Sphere 

fnL = I <5L'(f) ■ nr)d\ = 0, (24) 

J Svhere 



with 

«- = -XF r pr'iAn2ir) + 3Bn2ir))dr. (25) 
M Jo 

It is clear from Eq. ()24|) that toroidal modes do not enter into play, so we will work from now 
on only with spheroidal modes and we will drop the 5* for simplicity. 

Once the Riemann tensor and the h"{t) are known for a gravitational wave coming from 
a generic direction {9, 0), it is easy to compute Qnia'i't) using Eq. lfT^ . The solution of Eq.((TJ 
can be written: 

oo oo 

u{f, t) = Y: un{r, = E — E <^n2.(r)g^2(i)- (26) 

ra=l n=l ^"2 a 

For spheroidal modes, Bi\i{t) = Bn,i,a{t), as defined in the previous paragraph. In this case, 
and its equation of motion (fT3j) is 

= 5i^2anha{t)- (27) 
The energy stored in the mode a at frequency Uni is (from Eq. (|THjl ): 

1 f+OO 

Es{nJ,a) = -I / Fr,Ut)e-'^"''dt\^ (28) 

^ J —oo 



Using Eqs.jini) and Eq.(j2HI) becomes 

1 /■ + 00 

E,{n,l,a) = -al6i,2\ / h^''\t)e~'''"''dt\'. (29) 

/ J —oo 

The only contributions have, as expected, 1 = 2. For these, we fix the value of n and replace 
/i" using Eq.(jni). 

For the a = mode, we have 

D J-OO V 5 
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thus 
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Similarly, for the other modes, we get: 
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+ 2 cos 6* sin cos 0(1 + cos^ ^^)(cos2 — sin^ 0) 3fJ(/;.+/;,* )]. 

We have assumed that h^^^{t) and h^^y^{t) decrease fast enough at ±00 to allow the 
integration by parts, h^{f) and hy.{f) stand for the Fourier transform of h^{t) and hx{t). 
In this paper we use the notation d^u) = dtait)e^^'^^ for the Fourier transform and 
a{t) = /I*"^ (iu;a(u;)e"'"*'^* for the inverse. In Eq. ljH^ all Fourier transforms are evaluated at 
the frequency /„2 = i^n2/2vr. 

If the source emits randomly polarized radiation, we can rewrite hx and h+ in terms of 
one possible intrinsic polarization of the source (Eq.(^ni)) and average over ip. We finally 
obtain 



(32) 
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As before, all Fourier transforms are evaluated at the frequency ujn2- 

It is important to stress at this point that if we sum the energies of Eq. we obtain a 
constant value: 

27r 

E,{n,2,tot) = —a'^u^^^{\~e+{u2n)?+\~eAi^2n)?)- (34) 

This means that the total sensitivity of our detector is always the same, independently of 

the GW direction. This is to be expected because of the symmetry of the detector, and it 

is one of the most important properties differentiating spherical detectors from cylindrical 

ones. Together with its greater sensitivity (keeping the same dimensions), this feature makes 

the new generation of resonant detectors much more advantageous than bars. 

Figure 1 shows, for each mode n, how the energy stored in each a varies as a function of 

9 and 0, the angles that give the direction of the incoming unpolarized GW. The rate 

15E,(n,2,a) 
7ran^n2(|e+p + |exP) 

is zero where the graph is blue, and maximal (|) where the graph is red. 

The total sensitivity is the sum of contributions coming from the five modes a and the 

sensitivity of each mode depends on the orientation of the source with respect to the chosen 

frame. This means that for different angular positions of the source with respect to the 

detector reference frame, there will be different distributions among the five modes of the 

total energy transfered from the GW to the sphere. In contrast, if the source is fixed in 

space, and the sphere moves because of the earth's rotation, the energy ratio stored in each 

of the five modes will be a periodic function of the arrival time, with period one sidereal 

day, while the sum of the five contributions will be time-independent. 
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The expected sources for GW resonant detectors are traditionally distinguished into: 1) 
burst of relatively great intensity, given by cataclysmic events as binary cohoaleshing, or su- 
pernovae explosions; 2) sources which emit GW bursts occasionall y b ut repeatedly in time 
and with peculiar statistic features, the so called " GW-bursters" 3) periodic sources, 
of longer duration and smaller spread in frequency, as pulsars or inspiraling binaries. 
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Such astrophysical objets are present in galaxies, and as the energy hitting a detector de- 
creases as the distance of the source squared |l7|, it is then reasonable to expect a GW 
signal coming from our galaxy to be more easily detected. 

Once the position of the source is known in galactic coordinates {b,lg), a set of space 
rotations gives the position of the source in the detector frame {xd,yd,Zd)- These rotations 
are derived and explained in detail in the appendix of ^|. For a given (6, Ig) one can then 
easily find the precise position {6{t),(f){t)) of the source with respect to the detector frame, 
with sidereal time t varying from to 24. 

As example, figure 2 displays the energy rates ea stored in the five modes a of a sphere near 
Leiden, Holland ( latitude / = 52.16°A^ , longitude L = 4.45°i?, as functions of Greenwich 
sidereal time (in sidereal hours) and of galactic longitude b for a randomly polarized source 
lying in the galactic plane {Ig = 0), where stars population is more important (and then 
their remnants population also). As before, red zones indicate maximal values. 
Finally, the example of a section (at 6 = 0, i.e. for a source near the GC) for the a = 
mode is represented in Figure 3. These graphics can be specially interesting for sources of 
type 2 and 3 mentioned above, for which a signal from the same location persists in time. 



IV. THE SPHERICAL GW DETECTOR WITH A SET OF RESONATORS ON 
THE SURFACE 

Until now we have looked at the spherical detector as a unique object interacting with 
GW. Nevertheless, because of the smallness of h, the magnitude displacement of the detector 
as a result of this interaction is tiny, of the order of hL, L being the detector size. 

For this reason small resonators of equal mass rrir « M and elastic constant kr are 
put on the surface of the GW resonant detector. The small resonators are tuned to the 
fundamental frequency of the detector, in order to amplify displacements of its surface at 
that chosen frequency, which for the sphere turns out to be 2. In this way, frequencies 
other than the fundamental mode and nearby frequencies will not be amplified by the system. 
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A. The general case of K resonators 



Imagine having K such resonators on the surface of the sphere. We are interested in 
what happens when a GW interacts with such a system, in the ideal case of high SNR (i.e. 
where noises are negligible). 

In line with previous work, in deriving the equations of motion for such a system we will 
use the formalism and notations of [5] , where the problem has already elegantly been posed, 
and a proposal for a six-detector configuration, the TIGA configuration, is presented. 
The equation of motion (Q) for small displacements now becomes 

P^^^ - P^Mr, t) - (A + /i) V( V ■ u(f , t))=Y: /(^) (r)^?^'') (t) (35) 

where 77 is a parameter that takes the values 0, Ic, Is, 2c, 2s, 1, 2, K = {a, i}. It is 
clear that the first components of this force density. 

Few = E r(^1^?"W = E rirlh'^iil (36) 

a a 

are the usual ones (see Eq.(jll)) attributable to GW, while 

Fres = j:fHr)9^Ht) (37) 



i=l 

are new contributions given by the elastic terms of resonators. 

Let us focus now on one resonator, which we assume to have only radial motion, linked 
with a tuned spring to the surface of the sphere, say at the position R^. We define the 
scalar-time dependent quantity 

Zi{t) = u{Ri,t) ■ h (38) 
as the radial projection of the surface displacement at this point. 

We also define qiit) to be the relative distance between the resonator and the sphere surface. 
Note that qi{t) is not an inertial coordinate, but Zi{t) + qi{t) is. Then the equation of motion 
for this last quantity is 

mriUt) + ^iit)) = -Kqi{t) + F^^{t), (39) 

where F^"^ is the radial component of the GW force density acting on the resonator, which 
will be neglected as we know that m^. << M. 
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Since 

g\t) = q,{t) , f{f) = kr6\f-Ri)Ri, (40) 
the force density component (j?fj) can be written: 

Fres = E krq^it)6^ir- R,)R,. (41) 

i 

Generalizing the same method used above (in particular, see Eqs. (fTT|) . (fT^ and (fT^ ), 
the solution of Eq. (jHS)) can be expressed as: 

V n,l,a 

= Bn,iAi)KiA^ (42) 

n,l,a 

with 

M J Sphere 

QZ,iJt) = t9'{t')sm{uj^,Ut-t'))dt'. (43) 
Jo 

The last equahty of Eq. defines 

Bn,lAt)=T.<iJll,aQll,a{t). (44) 
V 

In our case we get 

a'n,lSa,a' H T] = a' 

[ !§A^,i{R)Yn,iA^,,(l),) iir] = t. 

It is then straightforward to derive the equations of motion for the sphere modes Bn^iA^), 
by inserting Eq. lj^^ into Eq. ljHSj) : 



v 

K 



1=1 



(46) 



When noises are negligible, the motion is dominated by quadrupole modes (/ = 2, a = 
0, Ic, Is, 2c, 2s), the only ones interacting with GW. We limit ourselves to the n = 1 modes, 
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for which the cross-section is the most significant Note that for / = 2, a„^2 = o„ of 

Eq.dnnD. 

From now on, we simphfy the notation, writing 

5i,2,a(t) = B^{t) , ai,2 = a , Ai,2(i?) = A 

Yi,2A0i,(f).) = Y^{e,,<f),) = Y^,. (47) 

This last equahty defines the 5 x K matrix Y. 

Muhiplying for the sphere mass M, and caUing the five equal quantities Muf 2a = = 
kg, Eq. (j46|l can be rewritten 



MB^{t) + ksB^{t) = Mah'^{t) - ^ krAY^iqi{t). (48) 

i 

Defining now the vectors 

{qm} = q(t) , {Mah^it)} = = H(t) , 

{B^m = h{t) , (49) 

we can again simplify the expression for these equations: 

Mishit) + kslMt) - KAYciit) = H(t). (50) 

It is evident that b and h are 5-component vectors, q is a A'-component vector and I5 the 
identity matrix 5x5. Thus Eq. lj^UI) is a system of 5 equations. 

In the same way we can rewrite Eq. ()39j) in matrix form. To do that, we need to express Zi 
as a function of Ba- Using the approximation that only the n = 1 modes of the sphere are 
excited at high SNR, from Eq. ljnHj) we have: 



Zi{t) = u{Ri,t)-h= E B^it)AYa^. (51) 

n=l,l=2,a 

Defining the vector {zi}(t) = z(t), 

z{t) = AY^h{t). (52) 

Then Eq. flHTIjl becomes 

mr{q{t) + AY^h{t)) + Kqit) = 0. (53) 
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The K equations for q ()53p and the 5 equations for b ()50p can be rewritten as a system of 
K + b equations 
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(54) 



These are the equations of motion for a system composed of a homogeneous sphere with 
K resonators moving radially on its surface, all interacting with a GW. 

In the remainder of this paper, we will solve these equations once a configuration for the 
positions of the resonators is given. We first look at what happens in the case with only 
one resonator (A; = 1), in order to test the equations by looking at a simple example. Then 
we move on to the TIGA configuration {K = 6) proposed by Merkowitz and Johnson in P| . 



B. The case of one resonator at 9 = 

In this particular case we choose to put only one resonator on the surface of the sphere. 
Obviously each position is equivalent, but for simplicity let us put the resonator at the north 
pole of our coordinates system (or, equivalently, at the south pole because of the symmetry 
of the problem). 

In this case, the only nonvanishing component of the 5x1 matrix Y is Yqi = = y. 
Then Eq. (j54|l can be written: 
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(55) 

It is clear from Eq. (j55p that the only mode of the sphere which couples with our resonator 
is Bq, while all the others behave as if the resonator was not present. Their equations of 
motion are the same as given above in Eq. ()27|) . 

So we can reduce our system of equations to a system of 2 coupled equations with 2 un- 
knowns, taking only the first and the last lines of (j^ : 

Bo{t) + uoIBq = LolAyeq + a/i°(t) 

q{t)+^Uit) = -AyBo{t). (56) 

In Eq. fl56p we consider that the resonator frequency is perfectly tuned with the funda- 
mental frequency of the sphere, i.e. that 

M ~ m 

and we call e the (small) ratio between the resonator mass and the sphere mass. 
In order to solve Eq. (j3H|) we move to Fourier coordinates, giving: 



-uj a 



{uq — cj^)^ — Lu^ujQA'^y'^e 



auj'h\uj)- 



ijj — ijJn 



(57) 



As expected, the resonant frequencies ijJ± for both the coupled quantities are the same, 
and they are the splitting of the original frequency ujq. 
They can be found setting the denominators of Eq. ()57p equal to zero: 



UJ. 



± 



UJc^ 



\±Ay 



\ 



1 + 



A^yH 



+ 



(58) 
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At the frequency c^o we have Bo{tuo) = 0, and q{uJo) = -^h^{u!o) is proportional to 
the ratio between the sphere mass and the resonator mass, so we will have the greatest 
oscillations if we make the resonator mass as small as possible. 

Note that here we neglected the damping of both the coupled oscillators. This is because 
we are interested in burst signals, where the time-scale for damping an oscillation of the 
system is much longer than the time duration of a signal. Taking account of damping means 
changing Eq. ()57j) to cancel divergences which arise at u = uj± and make them smoother. In 
particular, Bq{ujq) will no longer vanish. Computations adding damping terms are analog 
and easily done. Figure 4 shows a plot of |g(co')p with a sphere damping of 10""^, a resonator 
damping of lO'^, uo = lO^Hz, e = 10"^ and h°{uj) = 10~^^Hz-\ 

The quantity we are really interested in is the one we can measure, q{uj), or its squared 
absolute value. Rewriting to show its dependence on and on the wave direction 

(see Eq.®) we get: 



\q(uj)\'^ = 7y sin^ 9\h+(uj) 

= f3{uj)sm^e\h+{u)\\ (59) 



or, assuming that the GW is unpolarized, and averaging over polarizations. 



/3(c.)(|e~+|2 + |e~xP)/2 

with 



sin^(^), (60) 



This turns out to be the same angular sensitivity obtained for the cylindrical bar , for 
which one finds that sensitivity is optimal for a plane orthogonal to the bar and drops with 
a factor sin^ 9 too. 

Having only one transducer on the sphere surface is a simple example showing that the way 
resonators are positioned can result in the loss of the main initial properties given by the 
spherical symmetry of the detector. In order to keep these properties as far as possible, we 
have to put more resonators on the sphere in appropriate positions. 
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C. The TIG A configuration 



With the hope that, in future, we will be able to place resonators moving radially on the 
sphere surface, and that we will be able to measure experimentally the displacements q of 
these masses, and assuming that all other interactions with the sphere-resonators system are 



negligible [18| with the exception ofgravitational waves, in this section we find the response 
of (i{uj) in the TIGA configuration 5| to a gravitational wave coming from a generic direction 
{6, 0). In order to do so, we have to solve the equations of motion ()54|) in Fourier space. 

In the TIGA configuration there are a total of six resonators: three of them have azimuth 
angle 6i = 62 = 63 = 6a and 0i = 0, 02 = 03 = ^ and the other three have 64 = 6^ = 
9q = Ob and 04 = |, 05 = vr, 06 = The angles 6a and Ob are between and 7r/2 and 
they are solutions of the equation 

45 cos^ ^ - 30 cos^ + 1 = 0. 

With this information, the Y matrix for this configuration is easily found: 



47r 



47r 
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3+v^ 
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3+v^ 
247r 



3+V5 
Sit 



3-VE 



3+v^ 







3-^5 
247r 



3-v^ 
8n 



3+V5 
2An 



3+V5 
8tt 



(62) 



It is important to recall a very usefull property of this matrix: YY^ = is proportional 
to the identity 5x5 matrix. Nevertheless, the inverse is not true: Y^Y has all diagonal 
values equal to ^ and all non diagonal values equal to It is easily seen that sums over 
lines or columns are always zero. Actually, Y-^Y is a 6 x 6 matrix obtained from two 5x6 
matrices, so it is clear that it cannot be invertible. 

Using this property it is straightforward to solve the equations for the so-called 
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mode — channels 



p{uj) = Yq(^) 



(63) 



and for h(uj) 



Pa[UJ) 



BJuj) 



2lT , 



The denominators vanish for the two spht frequencies 



1 ± — W— eA2 + —— + —eA^ 



Arc 



so one can rewrite Eg. ()fi4j) as 



Pa[UJ) 



Aauj^h''{uj) 



27r (u;2_^2)(^2_^2) 



with 



A. 



A_ 



—u^Ah'^, 
2n 



2 2 
UJl — UJl 



(64) 



(65) 



(66) 



(67) 



(68) 



UJ^ — ul 

We then have an expression giving a hnear combination of resonator displacements as 
a function of the spherical components of a gravitational wave, or we have each spherical 
component of /i as a function of different combinations of displacements. 

Nevertheless we are also interested in the inverse relation: we want to write each (ji as 
a combination of the h°'. In order to do that we have to express 6 quantities as a function 
of 5, which means inverting the matrix Y, a 5 x 6 matrix. We then need a supplementary 
condition for q. The condition to be added is true for the TIGA configuration only, in the 
case of a high SNR. In that case, the sum of all the six displacements is zero, and this can 
be seen by looking at the only non-trivial solution of Yq = 0. This means that we can add 
a sixth line to the Y matrix, a line made up of six unities, and a sixth null component to 
the p vector in Eq.(|63|). We will call the new 6x6 matrix Y and the new 6- vector p. 

Y is invertible, and we can easily find q using q = Y^^p, getting: 



37T 
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3 ' ' ' U 
+ ^(v/5-l)p2.H; 



+ y^{V5 + i)p2s{u;y, 

o 



96 (^) 



6 



(V5 + 1)P: 



(69) 



18 ^ ^-^^ ^ 6 ' . I 



Taking the squared absolute value, and rewriting h'^{uj) as a function of /i+ , and of 
the direction (^,0) (see Eq.(jni)) we can find the explicit behavior of the amplitude of the 
oscillations for each resonator as a function of the GW direction and polarization. 

Supposing the source to be randomly polarized, we can replace /i+ x with the expressions 
in Eq. lfTUI) . and average over polarizations (as we have already done for the sphere without 
resonators). This is equivalent to replacing 



(70) 



The ratios 



with [3{uj) 



Ai{e, $) 







MP 




/3M(|e 
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ex 


2)/2 



2_^,i"m,^,2_,.,2 \2 , give the sensitivities of the six resonators as a function of 



the direction of GW provenance. 

The explicit expression for one of the A.i previously defined is: 
'V3 . 2, 73 



A, 



V 3 1 

sin^ 6 H (v^ + l)(cos^ sin6' COS0) H — (a/S + 1) (cos ^ sin 6 sin< 

6 18 6 
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a/3 1 1^ 

--7^(V5- l)(l + cos2^)(cosV- 1/2) - -(VS- 1)(1 + cos^ 9) cos sin 



+ 



V 3 1 V 3 

(Vb — l)(sin6' sin0) H — [Vb + l)(sin 9 cos0) H (Vb — 1) cos cos sin ( 

18 6 9 



— -(v^ — 1) cos6'(cos^ — sin^ • 
6 



(71) 



The other five Ai have analogous expressions. Their angular dependence is shown in Figure 
|SJ while Figure 6 shows the angular dependence of their sum. 

Note that the sum of the six sensitivities is no longer a constant. This is because, 
as we noted before, putting resonators on the surface of the sphere means destroying the 
original spherical symmetry. Nevertheless, the total sensitivity is never zero, so for the 
direction of a GW from any provenance, there is at least one resonator which will be excited 
significantly. Omnidirectional sensitivity is then preserved in the TIGA configuration, even 
if omnidirectional symmetry is broken. 

In Figure U\ we take into account the earth's rotation, as explained in the appendix of 
for a fixed source near the galactic center. Angular sensitivities for the six resonators 
and for their sum are shown as a function of sidereal hour and galactic longitude. As before, 
we take as an example the sphere near Leiden, Holland ( latitude / = 52.16°A^ , longitude 
L = 4.45°E). 



V. SUMMARY AND CONCLUSIONS 



In this work, the general problem of the interaction of a spherical detector with GWs has 
been analyzed, in both the cases of a simple elastic sphere and of a system sphere + trans- 
ducers. The five responses of the quadrupole fundamental modes and of resonators motions 
for the TIGA configuration have been found, explicitly identifying angular dependence and, 
for GW radiation coming from a source in the galactic plane, sidereal time dependence as 
well. 

We find the expected omnidirectional constant sensitivity for the ideal case of a simple 
sphere. This omnidirectional sensitivity is shown to be preserved in the case of the sphere 
equipped with TIGA transducers, nevertheless it is interesting to note that it is no longer a 
constant. This confirms the fact that putting resonators on the surface of the sphere means 
choosing some favorite directions, breaking the omnidirectional initial symmetry. 
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For current experiments, the results of this paper allow one to predict exactly how the 
resonators should behave once the position of a candidate source is known. Vice versa, rela- 
tions between the resonators' amplitude and {6{t), (pit)) can be inversed and whenever there 
is a statistically significative reason for thinking that a gravitational wave signal has been 
received, one can estimate the source direction by looking at the ratio between resonators' 
excitations and at the GW arrival time. 

Nevertheless, all this work is based on an ideal sphere and ideal radial moving resonators. 
Deviations and errors on the above predictions exist in a real spherical detector experiment. 
In particular, the asymmetry caused by the hole drilled through the center for the suspension 
system induces a degeneracy of the quadrupole five fundamental modes, which in is 
estimated to be less than 1%. We have not derived the consequences of this degeneracy 
in this paper, but it would be interesting to do that in order to better match theoretical 
predictions with experiments. 



[1] www.minigrail.nl . 

[2] www.das.inpe.br/ graviton . 

[3] C. Z. Zhou and P. F. Michelson, Phys. Rev. D 51, 2517 (1994). 
[4] J. A. Lobo and M. A. Serrano, Europhys. Lett. 35, 253 (1996). 
[5] S. M. Merkowitz and W. W. Johnson, Phys. Rev. D 51, 2546 (1995). 
[6] S. M. Merkowitz and W. W. Johnson, Phys. Rev. D 56, 7513 (1997). 

[7] R. V. Wagoner and H. J. Paik, Gravitazione Sperimentale, Accademia Nazionale del Lincei, 
Rome, 1977. 

[8] N. Ashby and J. Dreitlein, Phys. Rev. D 12, 336 (1975). 

[9] L. Landau, E. Lifshitz, Theorie de I'elasticite, Moscou, 1960. 
[10] C. W. Misner, K. S. Thorne, J. A. Wheeler, Gravitation, W. H. Freeman & Company (1970). 
[11] J. A. Lobo What can we learn about gravitational wave physics with an elastic spherical an- 
tenna?, Phys. Rev. D 52, 592 (1995). 
[12] L. Landau, E. Lifchitz, vol.1 Mecanique, Moscou, 1960. 
[13] E. Coccia and J. A. Lobo and J. A. Ortega, Phys. Rev D 52, 3735 (1995). 



24 



[14] K. S. Thorne, in 300 Years of Gravitation, edited by S. Hawking and W. Israel (Cambridge 

University Press, Cambridge, England, 1987). 
[15] S. M. Merkowitz, Ph.D. thesis, Louisiana State University (1995). 

[16] F. Dubath, S. Foffa, M. A. Gasparini, M. Maggiore and R. Sturani Experimental signature of 

gravitational wave bursters Phys. Rev. D 71, 124003 (2005). 
[17] E. Coccia, F. Dubath, M. Maggiore Possible sources of gravitational waves bursts detectable 

today Phys. Rev. D 70, 084010 (2004). 
[18] In particular, we consider the case where another amplification device, such as SQID, is NOT 

tuned to the frequency of resonance of the sphere. 

Acknowledgments 

I would like to thank Florian Dubath and Stefano Foffa for many helpful discussions. 
A special thank to Angela Haden, Danielle Chevalier, Andreas Malaspinas, Marti Ruiz- 
Altaba and Fabio Dubath. 

This work is partially supported by the Swiss National Fund (FNS). 



25 



5 11 15 ^ , , 5 10 15 ;p 




e2s 

5 10 ' 15 20 




FIG. 2: Sensitivity of the five fundamental modes as a function of the GW source position in the 
galactic plane, (6), and of sidereal time in sidereal hours, (t), for a sphere near Leiden, NL (latitude 
/ = 52.16°A^, longitude L = 4:A5°E). As in Figure 1, red zones indicate maximal values (Color 
online) . 
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FIG. 3: eo of Figure 2 as a function of t in sidereal hours for b=0 (source near the galactic center) 




omega 

FIG. 4: with a sphere damping of 10~^, a resonator damping of 10^^, ujq = W^Hz, e = 10~ 

and h°{uj) = lO'^^Hz-K 
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FIG. 5: Sensitivity of the six resonators in TIG A configuration as a function of the GW arrival 
direction {9,(j)). As in Figure 1, the ratio Ai is zero where the graph is blue, and maximal where 
the graph is red (Color online). 
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FIG. 6: The same as Figure 5 for the sum of the Ai (Color onUne). 
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FIG. 7: The same as Figure 5, but as a function of sidereal hours (t) and galactic longitude (b) of 
the GW source, for a spherical detector near Leiden, NL, / = 52.16°A^, L = 4.45°-E (Color online). 
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FIG. 8: Sum of the Ai of Figure [7| as a function of sidereal time for b=0, i.e. a source near the 
galactic center (Color online). 
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